Signature of the interaction between dark energy and dark matter in galaxy clusters 



Elcio Abdalla^, L. Raul Abramo'*, Laerte Sodre Jr.*^, Bin Wang'^ 

"'Instituto de Fisica, Universidade de Sao Paulo, CP 66318, 05315-970, Sao Paulo, Brazil 
^ Departamento de Astronomia, Instituto de Astronomia, Geofisica e Ciencias Atmosfericas da Universidade de Sao Paulo 
Rua do Matdo 1226, Cidade Universitdria, 05508-090, Sao Paulo, Brazil 
'^Department of Physics, Fudan University, 200433 Shanghai, China 



Abstract 



o 
o 

(N 



We investigate the influence of an interaction between dark energy and dark matter upon the dynamics of galaxy clusters. 
We obtain the general Layser-Irvine equation in the presence of interactions, and find how, in that case, the virial theorem 
stands corrected. Using optical, X-ray and weak lensing data from 33 relaxed galaxy clusters, we put constraints on 
the strength of the coupling between the dark sectors. Available data suggests that this coupling is small but positive, 
indicating that dark energy might be decaying into dark matter. Systematic effects between the several mass estimates, 
however, should be better known, before definitive conclusions on the magnitude and significance of this coupling could 
be established. 
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1. Introduction 

Cosmological accelerated expansion is by now a well- 
established observational fact [U El [3] , leading either to an 
asymptotically de Sitter cosmology, plagued with an aston- 
ishingly small cosmological constant, or else to a universe 
filled up to 80% with a strange dynamical component with 
negative pressure - dark energy 0]. 

If dark energy contributes a significant fraction of the 
content of the Universe, it is natural, in the framework of 
field theory, to consider its interactions with the remaining 
fields of the Standard Model and well- motivated extensions 
thereof. For lack of evidence to the contrary, interactions 
of dark energy or dark matter with baryonic matter and 
radiation must be either inexistent or negligible. Never- 
theless, some level of interaction between dark energy and 
the dark matter sector, which is present in most extensions 
of the Standard Model, is still allowed by observations. 

The possibility that dark energy and dark matter can 
interact has been studied in [S]-[TS], among others. It 
has been shown that the coupling between a dark energy 
(or quintessence) field and the dark matter can provide a 
mechanism to alleviate the coincidence problem [51 [TD] . A 
suitable choice of the coupling, motivated by holographic 
arguments, can also lead to the crossing of the phantom 
barrier which separates models with equations of state 
w — p/p > —1 from models with w < —1 [11] - see also 
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[T51 [13]. In addition, it has been argued that an appro- 
priate interaction between dark energy and dark matter 
can influence the perturbation dynamics and affect the 
lowest multipoles of the CMB spectrum, accounting for 
the observed suppression of the quadrupole [TH [TS] . The 
strength of the coupling could be as large as the fine struc- 
ture constant [T5] [TB] . Recently, it was shown that such 
an interaction could be inferred from the expansion history 
of the universe, as manifested in, e.g., the supernova data 
together with CMB and large-scale structure [17]. Never- 
theless, the observational limits on the strength of such an 
interaction remain weak |18j . 

A complementary and fundamentally different way in 
which the coupling between dark energy and dark mat- 
ter can be checked against the observations is through its 
impact on large-scale structure. If dark energy is not a 
cosmological constant, it must fluctuate in space and in 
time - and, in particular, if dark energy couples to dark 
matter, then it must surely be dynamical. If that is the 
case, dark energy affects not only the expansion rate, but 
the process of structure formation as well, through den- 
sity fluctuations, both in the linear [TU], [IS]- [22] and the 
non- linear [23l [23] regimes. The growth of dark matter 
perturbations can in fact be enhanced due to the coupling 
between these two components [TU [T31 |2S] . 

Recently, it was suggested that the dynamical equilib- 
rium of collapsed structures would be affected by the cou- 
pling of dark energy to dark matter, in a way that could 
be observed in the galaxy cluster Abell A586 [25]. The 
basic idea is that the virial theorem is distorted by the 
non-conservation of mass caused by the coupling. 

In this paper wc show precisely how the Layser-Irvine 
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equation, which describes the flow to viriahzation [27j . is 
changed by the presence of the coupling, in such a way that 
the final state of equilibrium violates the usual virial con- 
dition, 2K+U — 0, where K and U are respectively the ki- 
netic and the potential energies of the matter constituents 
in an isolated system. We show that this violation leads 
to a systematic bias in the estimation of masses of clus- 
ters if the usual virial conditions are employed. Although 
it is still possible that systematic errors from observations 
smear the results, the fact that some shift in the mean 
value of the coupling for two independent sets of observa- 
tions (compared to the third set) may signalize some new 
physics. 

Even though the uncertainties associated with any in- 
dividual galaxy cluster are very large, by comparing the 
naive virial masses of a large sample of clusters with their 
masses estimated by X-ray and by weak lensing data, 
we may be able to constrain such a bias and to impose 
tighter limits on the strength of the coupling than has 
been achieved before. 



2. PHENOMENOLOGY OF COUPLED DARK 
ENERGY AND DARK MATTER MODELS 

Quite generically, at the level of the cosmological back- 
ground an interaction between dark matter and dark en- 
ergy manifests itself as a source term in the continuity 
equations of both fluids: 



Pde 



Pdrn + "iHpdm = V' 
"iHpdeil+Wde) = -V' , 



(1) 



where a dot denotes time derivative, H is the expansion 
rate, pdm and pde are respectively the energy densities of 
dark matter and dark energy, and Wdm and Wde are their 
equation of state parameter. Notice that the continuity 
equation still holds for the total energy density pTot — 

Pdm + Pde- 

Phenomenologically, one can describe the interaction be- 
tween the two fluids as an exchange of energy at a rate 
proportional to the total energy density [51 [TT]: 



^ = CHp 



'Tot 



(2) 



We are interested in collapsed structures - places where 
the local, inhomogeneous density a is far from the aver- 
age, homogeneous density p. In that case the continuity 
equation for dark matter reads: 

(^dm + 3H(7dm + V (cTdrnVdrn) = (H {(J dm + 0"de) , (3) 

where Vdm is the peculiar velocity of dark matter particles. 

In this work we will consider the local density of dark 
energy to be proportional to the local density of dark mat- 
ter, Ude = bem'^dm- If for a givcu modcl the dark energy 
component is very homogeneous, hem ~ 0. We do not con- 
sider the case where h^m. depends on the size and mass 



of the collapsed structure - although this should proba- 
bly happen in realistic models of structure formation with 
dark energy [24' . Hence, the continuity equation with dark 
matter coupled to dark energy reads: 



O'dm + SHadm + V {cfdmVdm) — (Hudn 



(4) 



where C = C(l + ^em) is the effective coupling in a virialized 
structure. Notice that different dark energy models pre- 
dict different levels of dark energy perturbations |24l |25] , 
hence any constraints we derive from observations of col- 
lapsed structures will be in some sense degenerate with the 
perturbative properties of the dark energy sector. 



3. LAYZER-IRVINE EQUATION IN THE PRES- 
ENCE OF COUPLING 

We will use Newtonian mechanics to derive equilibrium 
conditions for a collapsed structure in an expanding Uni- 
verse. The acceleration due to the gravitational force is 
given by: 

(a Vdm)' = -a , (5) 

where a is the scale factor and tp is the (Newtonian) gravi- 
tational potential. Multiplying both sides of this equation 
by (JdmCiVdm, integrating over the volume and using the 
continuity Eq. (|4|, we get that the left-hand side becomes: 



(a^Kdm) ' - a^CHKdm 



(6) 



where the kinetic energy of dark matter is given by: 

Kdm ^\ J '"dm'^dmdV . (7) 

The right-hand side of the equation, on the other hand, 
becomes: 



(1 + Km) [-a' (Udm + HUdm) + KHa^Udr, 



(8) 



where we have used the Poisson equation, the fact that 
UTot = (1 + bem)o'dm, and the definition of the potential 
energy of a distribution of dark matter particles: 



Udn 



1 



:G 



{x)(Jdm{x') 



dVdV' 



\X — X' 



(9) 



The identity between Eqs. ^ and (|8| is the generalizaton 
of the Layzer-Irvine equation [27. describing how a col- 
lapsing system reaches a state of dynamical equilibrium in 
an expanding universe. One can see that the presence of 
the coupling between dark energy and dark matter changes 
both the time required by the system to reach equilibrium, 
and the equilibrium configuration itself. For a system in 
equilibrium {Kdm = Udm = 0) we get the condition: 



(2 - QKdm + (1 + bem){l - 2C)C/*,„ - 



(10) 



Taking C = b^m = we recover the usual virial condition. 
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4. MASS ESTIMATION AND LIMITS ON THE 
COUPLING 

Galaxy clusters are the largest virialized structures in 
the Universe, and their mass content is supposed to be rep- 
resentative of the universe as a whole - see, e.g., [2H]. They 
are composed of hundreds of galaxies, with the largest frac- 
tion of their baryonic mass in the form of hot. X-ray emit- 
ting gas - not stars. Clusters are widely believed to be to- 
tally dominated by dark matter [29j, and are conspicuous: 
existing surveys have already detected many thousands of 
clusters, and upcoming surveys will map much more. 

Cluster masses can be estimated in a variety of ways. 
Weak lensing methods use the distortion in the pattern 
of images behind the cluster (which acts as a lens) to 
compute the projected gravitational potential due to that 
cluster. Knowing the distances to the cluster and to the 
background images, one can derive the mass that causes 
that potential. An independent mass estimation can be 
obtained from X-ray observations if we assume that the 
ionized gas is in hydrostatic equilibrium. In this case the 
cluster mass can be determined by the condition that the 
gravitational attraction is supported by the gas pressure. 
Finally, wc can measure radial velocities and the projected 
distribution of cluster galaxies and, by assuming that clus- 
ters are virialized, one can infer their masses using the fact 
that U (X a"^ but K (x a, hence U/K (x M. 
However, Eq. 
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tells us that when there is coupling 
between dark matter and dark energy, the equilibrium con- 
dition depends on the coupling as: 



(1 + bern) 



Ud-n 

k7 



dm 



, l-C/2 
■ 1-2C 



(11) 



Hence, the mass that is estimated under the assumption 
that 1^ = is biased with the respect to the actual mass 
by a factor of (1 - C/2)/(l - 2C). 

One can compare directly the mass obtained through 
the virial theorem with that determined by other meth- 
ods. Notice that the total mass of a cluster is, within 
our approximations, the integral of ctto* = (1 + bem)o'dm 
over the volume; hence, the mass estimated through X- 
rays and weak lensing is related to the virial mass as 

M = {l + b,,„,)Md„r. 

Therefore, comparing the masses estimated through the 
conventional virial hypothesis with the masses estimated 
either by weak lensing or by X-ray data, we get: 



Mx = MwL = My„ X 



l-C/2 
1-2C 



(12) 



Hence, there are three tests one can make with these 
datasets: 



h = Mx/MwL ■ 



(13) 
(14) 
(15) 



If our interaction model is right, the first two tests, /i and 
/2, should agree with each other and put similar limits on 
the effective coupling parameter C, while the third test, 
/a, should only be a check of our method, and its value 
should be equal to one unless there are unknown system- 
atics affecting our mass estimates. Notice that either a 
violation of the equivalence principle for dark matter or 
a self-interaction of dark matter with itself, such as sug- 
gested by [30], could also be tested by comparing the dif- 
ferent mass estimates - see also [26] . 

To what extent the data currently available on cluster 
masses allows us to constrain the coupling parameter? To 
compare masses obtained with the different methods, we 
have assumed that the mass profile of the clusters is de- 
scribed by a Singular Isothermal Sphere (SIS). The main 
reasons are that weak-lcnsing mass estimations need to 
adopt a parametric model for the mass distribution in or- 
der to avoid the so-called mass-sheet degeneracy [HTj; this 
model is largely adopted in weak lensing studies, and most 
weak-lensing and X-ray mass estimations are possible only 
for radii significantly smaller than the virial. The main ad- 
vantage of this model is that it has a single parameter - 
the velocity dispersion along the linc-of sight a-y - which 
can be easily determined: directly from the observed ra- 
dial velocities in the virial estimation; from the X-ray tem- 
perature, [ax — kTx / {^-niH), where ^=0.61 is the mean 
molecular weight]; and from the fitting of the shear field 
in the case of weak-lensing. Since in this model the mass 
inside a given (projected) radius R is M{< R) = ira'^R/G, 
to compare the masses obtained by each method we need 
only to compare the velocity dispersions. 

For this exercise, we have analyzed data from galaxy 
clusters studied in [STJ 133] . Our sample has 33 clusters 
and was selected due to the homogeneity in the analy- 
sis procedure and avoiding clusters with evidence of dy- 
namical activity, like substructures. Ref. f31] presents a 
weak-lensing analysis of 24 galaxy clusters also observed 
in X-rays, verifying that clusters whith a hot inter galactic 
medium {Tx > 8 keV) are very active. For our analysis, we 
selected from this paper 14 clusters (from the 15 clusters 
with X-ray temperatures lower than 8 keV, Abell 1651 also 
has evidence of significant substructure [34]). The cluster 
A586, discussed in [32 , also seems to be in equilibrium. 
The remaining 18 clusters of our sample comes from [33]. 

We have used this dataset to test the theory that the 
usual virial mass is biased by a factor (1 — 2i^)/(l — (^/2) 
when compared to other mass estimates. Although the 
three datasets have asymmetric errors, we have assumed 
that the likelihood function associated with the three tests 
is symmetric, with width a = ^a^a^. With these as- 
sumptions, the likelihood function of test i is: 



Ni 

Li oc JJ^ exp ■ 

n=l 



1 



1^20 
1 - G/2 



(16) 



where the product runs over the data for each galaxy clus- 
ter n. If our model is correct we should get Ci = C2 = Ci 
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and (3 = 0. 

In the three panels of Fig. 1 we show the probabihty 
distribution functions for /i, /2 and /a, all computed for 
the top-hat prior —0.2 < Q < 0.2, which is more than suf- 
ficient to include the 3-cr limits —0.12 < C,i < 0.06 found 
by [18'. The shadowed regions in the top and middle pan- 
els mark the 67% and 95% Confidence Level (C.L.) limits 
for Ci and ^2- 

For both the /i and /2 tests we get a best-fit value of 
C ~ 0.03—0.04, while for the check /s we indeed get that Cs 
is consistent with zero with a high statistical significance. 
Our 95% C.L. hmits are 0.0 < C < 0.06 for the test A and 
0-0 ^ C ^ 0-09 for /2. This compares favourably with the 
95% C.L. constraint -0.095 < C < 0.035 obtained by [IH]. 

For /i and /2 we obtain that the null hypothesis = 0) 
is marginally consistent with the data, at the edge of the 
95% C.L. region. The statistical improvement between 
the null hypothesis and the best-fit model with coupling is 
weak, though: we get a Ax^/d.o.f. ~ 0.2 for both the /i 
and /2 tests. 

We have also computed the Bayesian Information Cri- 
terion (BIC) [35, to weigh if, and by how much, a coupling 
is necessary. For the test fi we get Agj^ « —0.5, while for 
the test /2 we get A|j(^ « —2.1. We can also estimate the 
level of systematic uncertainties that would turn the BIC 
against our model (i.e., when Abic = 0): an enhancement 
of 20% of all uncertainties would make Agj^ « 0, but still 
^Bic ~ in order to make A|jq « it would take an 

enhancement of 70% of the uncertainties. 

The reliability of these constraints, however, are dis- 
putable, due to possible systematic effects in the mass de- 
terminations. For example, a virial mass estimate is af- 
fected by the assumptions about the galaxy orbits, cluster 
morphology, mass distribution, identification of interlop- 
ers, etc. [331 EHl [37] , and the robustness of our simple SIS 
model does not mean that it is insensitive to (unknown) 
systematics. We can have a hint on the impact of these ef- 
fects for the constraints on the effective coupling constant 
by increasing its possible range of variation. For this ex- 
ercise we have just redone the analysis assuming that the 
actual errors are twice the internal errors. In this case the 
original results (Ci = 0.029 ± 0.015 and (2 = 0.044 ± 0.027 
at the 68% C.L.) change to Ci = 0.029 ± 0.030 and 
C2 — 0.044 ± 0.037, i.e., the most probable value of Ci^2 
does not change (as expected in this case), but the error 
in the estimates almost doubles, reducing to about one 
sigma the level of detection of a non-zero coupling con- 
stant. Moreover, since we can only constrain the effective 
coupling parameter = (1 -I- &em)Ci the bias between dark 
matter and dark energy in virialized structures could en- 
hance (if bem > 0) or suppress (6em < 0) our ability to 
constrain the true coupling ^. 

These results show that the reliability of a detection of 
the effective coupling parameter requires very good knowl- 
edge of possible systematic errors. Nevertheless, it also 
shows that if in the future we can produce a sample with 
reliable mass estimates and controlled systematics, we will 



indeed be able to constrain 2 and verify whether this 
hint of a coupling between dark matter and dark energy 
found with current data is confirmed. 

5. CONCLUSIONS 

We have estimated the effective coupling between dark 
energy and dark matter through the internal dynamics 
of galaxy clusters. In the presence of coupling, the flow 
of mass and energy between the components changes the 
virial condition in a way that can be tested by comparing 
different estimators for the mass of clusters. We searched 
for this signature in 33 galaxy clusters for which reliable 
X-ray, weak lensing and optical data were available. 

Our results indicate a weak preference for a small but 
positive effective coupling constant ^ - in line with pre- 
dictions made by some of us [TTJ [T^ . Since the statistical 
significance is still low (Ax^/d.o.f. ~ 0.2), it is paramount 
that more clusters (with homogeneous mass determina- 
tions and good control of systematics) be tested. If a 
significant indication of such coupling is still found, this 
would open a tantalizing new window on the nature of the 
dark sector. 
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Figure 1: Normalized probability distribution functions for the tests 
/l (top panel), /2 (middle panel) and (bottom panel.) The shad- 
owed regions on the top and middle panels indicate the 67% and 95% 
C.L. limits. 
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